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^/s^ ^ 20 - 200 GeV: data versus parton energy loss models 
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Abstract. Experimental results on high transverse momentum (leading) hadron spectra in nucleus- nucleus 
collisions in the range y?s NN ~ 20 - 200 GeV are reviewed with an emphasis on the observed suppression 
compared to free space production in proton-proton collisions at the corresponding center-of-mass energies. 
The transverse-momentum and collision-energy (but seemingly not the in-medium path length) dependence 
of the experimental suppression factors measured in central collisions is consistent with the expectations 
of final-state non-Abelian parton energy loss in a dense QCD medium. 

PACS. 12.38.Mh - 13.87.Fh - 24.85. +p - 25.75.-q 



1 Introduction 

High-energy nucleus-nucleus collisions offer the only ex- 
perimental means known so far to concentrate a signif- 
icant amount of energy (O(TeV) at the RHIC collider) 
in a "large" volume, O(10 3 fm 3 ), under laboratory condi- 
tions. The expectations, based on first-principles lattice- 
QCD calculations |T], are that for energy densities above 
e w 0.7 GeV/fm 3 , hadronic matter will undergo a phase 
transition towards an extended volume of deconfined and 
massless quarks, and gluons: the Quark Gluon Plasma 
(QGP). The scrutiny of this new state of matter aspires 
to shed light on some of the open key questions of the 
strong interaction (confinement, chiral symmetry break- 
ing, structure of the QCD vacuum, hadronization) that 
still evade a thorough theoretical description due to 
their highly non-perturbative nature. 

The production of an extremely hot and dense partonic 
system in relativistic heavy-ion reactions should manifest 
itself in a variety of experimental signatures. One of the 
first proposed "smoking guns" of QGP formation was "jet 
quenching" 0. Namely, the disappearance of the colli- 
mated spray of hadrons resulting from the fragmentation 
of a hard scattered parton due to the "absorption" of the 
parent quark or gluon as it traverses the dense strongly 
interacting medium produced in the reaction. Extensive 
theoretical work on high-energy parton propagation in a 
QCD medium 4,5,6,7. has shown that the main mech- 
anism of parton attenuation is of radiative nature: the 
traversing parton loses energy mainly by multiple gluon 
emission ("gluonstrahlung"). This medium- induced non- 
Abelian energy loss results in several observable experi- 
mental consequences: 



(i) depleted production of high px (leading) hadrons 0] , 

(ii) unbalanced back-to-back di-jet azimuthal correlations [HI 

El, 

(iii) modified parton fragmentation functions (energy flow 
and particle multiplicity within the final jets) (1011111 

D2|. 

The most simple empirically testable (and most easily the- 
oretically computable) consequence of jet quenching is the 
suppression of inclusive high px hadron spectra relative to 
their production in proton-proton collisions in free space. 
Since most of the energy of the fragmenting parton goes 
into a single leading hadron carrying a large fraction of 
the original parton energy ((z) = Phadron/p P arton ~ 0.5 
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Abelian energy loss should result in a significantly sup- 
pressed production of high px hadrons 0] . The amount of 
suppression is proportional to two physical properties of 
the medium |5lfelT7) : 

(i) the initial parton (gluon) density, dN 9 /dy, or, equiv- 
alently, the transport coefficient q = (k T )/X (which 
measures the average transverse momentum squared 
transferred to the projectile parton per unit path length), 

(ii) the square of the traversed path- length, I? . 

In this contribution, experimental results on inclusive sin- 
gle hadron production at high pt in nucleus- nucleus (A+ A) 
collisions at top CERN-SPS energies (y/s w 20 GeV), in- 
termediate RHIC energies (y/s = 62.4 GeV), and maxi- 
mum RHIC energies (y/s = 200 GeV), are reviewed and 
confronted to the theoretical predictions (i) and (ii) of 
non-Abelian parton energy loss in a dense QCD medium. 
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In the absence of any initial- or final-state medium ef- 
fect, the total hard interaction probability in a given A+B 
collision is only due to independent parton collisions which 
add incoherently. Based on individual point-like scatter- 
ing and general QCD factorization arguments E2|, the 
total hard cross-sections in A+B collisions can be written 
as the incoherent sum of all possible interactions between 
partons of nucleus A and partons of nucleus B. Namely, 



Eda^i hx /d 3 p = A-B -Eda^ r J hx /d 3 p, 



(1) 



where A (B) is the mass number (i.e. the number of nucle- 
ons) of nucleus A (B). Direct experimental measurements 
of hard processes in nuclear collisions, such as Drell-Yan 
production in Pb+Pb at CERN-SPS Q3|, or prompt-7 Q21 
HB| and total charm yields JTJ in Au+Au at RHIC, sup- 
port such a scaling. For a given centrality bin, Eq. l(T)l 
translates into the so-called "N co u (binary) scaling" be- 
tween hard p+p cross-sections and A+B yields: 



/d 3 p(b) = (T AB (b)) 



Eda h p ° r J hx /d 3 p, 



(2) 



where Tab (b) (Glauber nuclear overlap function) gives the 
number of nucleon-nucleon (NN) collisions in the A+B 
transverse overlap area at impact parameter b. The stan- 
dard method to quantify the effects of the medium in a 
given hard probe produced in a A+A reaction is given by 
the nuclear modification factor. 



RAA(PT,y,b) 



"hot /dense QCD medium" 
"QCD vacuum" 
d 2 N AA /dydp T 

(T AA (b)) ■ d 2 a vv jdydp T '' 
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7i° spectra 0-10% central A+A: 



• PHENIX Preliminary AuAuNs = 200 GeV 
PHENIX AuAu,\ls= 130 GeV 
PHENIX Preliminary AuAu, \Js = 62.4 GeV 
▲ WA98 PbPb,\ls = 17.3 GeV 
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Fig. 1. High pt neutral pion spectra measured in central A+A 
collisions at different center-of-mass energies: y/s NN = 17.3 
GeV EH], 62.4 GeV HJJ, 130 GeV and 200 GeV |2^liM| . 



7T°- ± + 



which measures the deviation of A+A at b from an inco- 
herent superposition of NN collisions, at transverse mo- 
mentum px and rapidity y. 



2 A+A collisions at v /s^7 « 20 GeV. 

Three nucleus-nucleus experiments at CERN-SPS mea- 
sured hadron production above pt = 2 GeV/c. WA98 and 
CERES/NA45 measured tt° and tt* in Pb+Pb 18 and 
Pb+Au EH] reactions at y/s NN = 17.3 GeV resp., whereas 
WA80 measured tt° in S+Au at y/s^ = 19.4 GeV jUJ. 
At these relatively low center-of-mass energies, the cross- 
sections for hard-scattering are extremely low and the 
maximum transverse momenta measured was pt ~ 4 GeV/c 
in the three cases 1 . Nonetheless, the power-law tail charac- 
teristic of elementary parton-parton interactions is clearly 
apparent in the measured A+A spectra above pr ~ 2 
GeV/c (Fig.^l. Unfortunately, no baseline high pt pions 
were measured in p+p collisions at SPS at the same cm. 
energies as heavy-ions, and extrapolations from higher-^/s 
data were used to obtain the expected d 2 a pp /dpTdy spec- 
trum needed to compute the nuclear modification factor, 
via Eq. ©. 

1 Which, yet, roughly corresponded to a remarkable ~l/2 of 
the kinematical limit, p% ax = yfs/2 « 9 GeV/c. 



In [23 ! we compared several proposed p+p 
X parametrizations, to the existing data in the range 
y/s ~ 16 - 20 GeV, and found that the parametrization 
of Blattnig et. al 26; reproduced reasonably well, within 
~25%, the shape and magnitude of the experimental pion 
differential cross-sections below px « 4 GeV/c. Using this 
p+p reference, we obtained the nuclear modification fac- 
tors for central A+A collisions at SPS shown in Fig. [5] 
Hadron production below pt ~ 1 GeV/c falls, as expected, 
below R AA — 1 (the assumption of independent point-like 
scattering does not hold for soft processes), but high-pr 
hadroproduction is, within errors, consistent with scaling 
with the number of NN collisions. Such a result is at vari- 
ance with the factor of ~ 2 enhancement observed in high 
Pt pion production in peripheral Pb+Pb reactions at the 
same energies [T%ll2"5] attributed to initial-state pr broad- 
ening as observed in fixed-target p + A reactions ( "Cronin 
effect") |2I]- This result points to the existence of an at- 
tenuating mechanism that reduces the underlying Cronin 
enhancement from R AA > 2 down to values consistent 
with R AA w 1 in central A+A. Indeed, theoretical pre- 
dictions of high pt 7t° production in central Pb+Pb at 
SPS including Cronin broadening, nuclear-modified par- 
ton distribution functions (PDF), and final-state partonic 
energy loss in an expanding system with initial effective 2 
gluon densities dN 9 jdy = 400 - 600 [2E] reproduce well 
the observed nuclear modification factor (yellow band in 

Fig. m 

The interesting conclusion that a moderate amount of 
jet quenching is already present in the most central heavy- 



2 Note that at SPS energies, the initial medium is probably 
more "quarkonic" than "gluonic". 
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ion reactions at SPS would require, however, a direct (and 
accurate) measurement of the high px P±P pion spectrum 
at a/s = 17.3 GeV. Unfortunately, the default minimum 
collision energy at RHIC in the proton-proton mode is 
y/s f» 48 GeV [55], making it difficult to directly com- 
pare at RHIC high pr hadroproduction in A+A and p+p 
collisions at center-of-mass energies comparable to SPS. 
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■ Pb+Pb -> 7i°+X 7% central [WA98] f-Js = 1 7.3 GeV) 
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Fig. 2. Nuclear modification factors for pions produced at 
CERN-SPS in central Pb+Pb [TS], Pb+Au [T5], and S+Au [H] 
at y/s NN « 20 GeV obtained using the p+p parametrization 
proposed in |25j . compared to a theoretical prediction [28] of 
final-state parton energy loss in a system with initial gluon den- 
sities dN 3 /dy = 400 - 600. The shaded band at R AA = 1 rep- 
resents the overall fractional uncertainty of the data (CERES 
data |19| have an extra uncertainty of ±15% not shown). 



3 A+A collisions at y/s^ = 62.4 GeV. 

The study of the excitation function of high px hadron 
suppression between top SPS and top RHIC energies was 
the main motivation behind the dedicated Au+Au run at 
RHIC intermediate energies (y/s NN — 62.4 GeV) carried 
out in April 2004. PHENIX measured neutral pions in the 
range p T = 1-7 GeV/c [2T| PHOBOS [2Q| and STAR [3T] 
measured inclusive charged hadrons up to pr ~ 4.5 GeV/c 
and 12 GeV/c respectively. However, as in the SPS case, 
no concurrent p+p reference measurement was performed 
at yfs = 62.4 GeV, and the corresponding Au+Au nu- 
clear modifications factors were determined using p+p — > 
h , 7T° + X differential cross-sections measured at the top 
CERN-ISR energies (y/s = 62 - 63 GeV) in the 70s and 
80s. As discussed in [20], the existing large inconsistencies 
(up to a factor of ~3) among the different ISR ir° data 
sets can be greatly reduced by removing the direct-y and 
rj contaminations not subtracted from the original "un- 
resolved" 7T° measurements. By doing so, one can obtain 
an averaged d 2 a pp ^ 7r o /dpT<iy reference spectrum at y/s = 
62.4 GeV with uncertainties ±25% for the Raa denomi- 
nator. 



Figure shows the PHENIX and STAR preliminary 
nuclear modification factors for high pt tt [25 and [HJ 
in central Au+Au collisions at y/s NN = 62.4 GeV obtained 
using the p+p references discussed in |29U31| . Above pr ~ 
5 GeV/c, there are ^3 times less produced pions and in- 
clusive hadrons than expected from point-like scaling of 
the p+p cross-sections. The magnitude of the suppression 
can be reproduced by models that include parton energy 
loss in a system with initial gluon densities dN 9 /dy = 
650 - 800 j3"2Tl3'3TEi4l (GLV formalism 0, green band) or 
transport coefficients (q) « 7 GeV 2 /fm [3"5ll5rJ] (Salgado- 
Wiedemann quenching weig hts 0, yellow band). The range 
Pt ~ 1 - 5 GeV/c shows, however, a significant rise and 
fall of the Au+Au spectra compared to the scaled p+p 
reference 3 . Although part of this effect can be certainly at- 
tributed to the expected collective radial flow and Cronin [35] - 
recombination enhancement 3 at 62.4 GeV, there are 
also py-dependent uncertainties related to the relatively 
poorly known shape of the ISR p+p spectra in this inter- 
mediate pr range [55] . Clearly, a dedicated RHIC proton- 
proton run at this collision energy would help to reduce 
these uncertainties and better constraint the theoretical 
predictions of parton energy loss models. 
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Fig. 3. Preliminary PHENIX and STAR nuclear modification 

for 7T° 121) and h ± 131) obtained in central 

.0 tit- 



factors, Raa(pt 

Au+Au at y/s^ = 62.4 G eV using the p+p — > tt u , + X ref- 
erences discussed in [29131) . The data are compared to two the- 
oretical predictions for parton energy loss in a dense medium 
with initial gluon density dN 9 /dy — 650 - 800 |32| or, equiv- 
alently, transport coefficient (q) — 7 GeV 2 /fm |35I36) . The 
error "bands" around each data point indicate the ~25% sys- 
tematic uncertainty of the ISR p+p baseline spectra at y/s = 
62.4 GeV. The shaded band at Raa = 1 represents the overall 
fractional uncertainty of the data (absolute normalization of 
Au+Au spectra and (Taa) uncertainties). 



3 An effect which is larger for the charged hadrons than for 
7T° due to the observed enhanced (anti)proton production |37) . 
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4 A+A collisions at = 200 GeV. 

Undoubtedly, one of the most significant results from the 
first 4 years of operation at RHIC is the large high px 
hadron suppression observed in central Au+Au reactions 
at ^/s~^ = 200 GeV. Above pt ~ 5 GeV/c, pions [2U, eta 
mesons PU, and inclusive hadrons (h ) |40II41| show all a 
"universal" factor of ~5 suppression compared to the cor- 
responding T^yi-scaled proton-proton yields. Such a deficit 
is not observed for direct photons 15, 16 (Fig.^J). The ob- 
served hadron suppression remains constant as a function 
of pt up to the highest transverse momenta measured so 
far (px w 14 GeV/c for 7r° [2H, see Fig. 0} in agreement 
with the parton energy loss model predictions. In the first 
theoretical predictions 5 , approximations of the underly- 
ing Landau-Pomeranchuk-Migdal (LPM) interference ef- 
fect in gluon bremsstrahlung resulted in a logarithmic de- 
pendence of the quenching factor on the parton energy 
and, therefore, in a Raa(jPt) that (slowly) increased for 
increasing hadron p^'s in apparent contradiction with the 
data. However, (i) the use of a realistic energy distribution 
of the emitted gluons (rather than the mean value) |42| . 
(ii) finite kinematic constraints (in the energy loss and 
in- medium path length), and/or depleted nuclear PDFs 

(in the EMC region for Bjorken x w 2p ^) > 0.2 val- 

ues corresponding to px > 12 GeV /c) 28 and, (hi) the 
increasing power-law (local) exponent of the parton spec- 
tra with pt , all explain the effectively constant Raa 
evolution as a function of transverse momentum. 
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Au+Au (central collisions): 

■ Direct y (PHENIX) 
Inclusive h* (STAR) 

• ti° (PHENIX Preliminary) 

O PHENIX Preliminary r\ 
GLV parton energy loss (dN^dy = 1 100) 
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Fig. 4. Raa(pt) measured in central Au+Au at 200 GeV for: 
direct photons [15I16| . inclusive charged hadrons [40| . n° I24| . 
and n |24) compared to theoretical predictions for parton en- 
ergy loss in a dense medium with dN 9 /dy — 1100 s 2HJ. The 
shaded band at Raa = 1 represents the overall fractional un- 
certainty of the data (absolute normalization of spectra and 
{Taa) uncertainties). The baseline p+p reference of the direct 
7 Au+Au data is a NLO calculation whose uncertainties are 
indicated by the dotted lines around the points |16| . 



A robust prediction of non-Abelian parton energy loss 
calculations is the expected cx L 2 dependence of the av- 
erage energy loss as a function of the in-medium path 
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Fig. 5. Preliminary PHENIX nuclear modification factor, 
Raa(<P), for 7T° production above pr = 4 GeV/c as a function 
of the azimuthal angle (f> with respect to the reaction plane in 
5 centrality classes of Au+Au at ^Js NN = 200 GeV |43) : com- 
pared to parton energy loss calculations [35) for an azimuthally 
asymmetric system with average transport coefficient (q) m 14 
GeV 2 /fm (yellow band, encompassing the limits of two dif- 
ferent prescriptions for the sampling of the energy loss). The 
lines around the experimental data show the uncertainties in 
the reaction plane and Raa determination. 



length 0. Such a behaviour predicted for a static QCD 
medium, turns into an effective oc L-dependence in an ex- 
panding QGP 6 . An interesting way to experimentally 
test the L dependence of the energy loss is by exploiting 
the spatial azimuthal asymmetry of the system produced 
in non-central nuclear collisions. Indeed, due to the char- 
acteristic almond-like shape of the overlapping matter pro- 
duced in A+A reactions with finite impact parameter, par- 
tons traversing the produced medium along the direction 
perpendicular to the reaction plane ( "out-of-plane" ) will 
comparatively go through more matter than those going- 
parallel to it ("in-plane"), and therefore are expected to 
lose more energy. In general, the total path length along a 
given azimuthal angle (f> with respect to the reaction plane 
is L{4>) ss 1— (e/2) cos(2(f>) 03], where e is the eccentricity 
of the system. By looking at the suppression pattern along 
different <f> trajectories one can test the L dependence of 
the energy loss. PHENIX 031 has recently presented nu- 
clear modification factors for high pt 7r° in Au+Au colli- 
sions at y/s~^ — 200 GeV binned in <fr angle with respect 
to the reaction plane (determined with the Beam-Beam- 
Counters at high rapidities) . The resulting Raa (</>) curves 
(Fig. [SJ show clearly a factor of ^2 more suppression out- 
of-plane ((f> = n/2) than in-plane (<j> — 0) for all the central- 
izes (eccentricities) considered. Theoretical calculations of 
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parton energy loss (based in the quenching weights for- 
malism 0) in an azimuthally asymmetric medium |35| 
predict a significantly smaller difference between the sup- 
pression patterns for partons emitted at <p = and (f> = ir/2 
(yellow bands in Fig. [SJ. The discrepancy model-data is 
stronger for more peripheral centralities (with correspond- 
ingly larger eccentricities) and challenges the underlying 
in- medium path- length dependence of non- Abelian parton 
energy loss (a detailed discussion of this observation can 
be found in and/or points out the necessity of an ad- 
ditional source of azimuthal anisotropy in pion production 
at high pt- Likely, collective elliptic flow is responsible for 
the extra boost of in-plane pions (even at pr values above 
4 GeV/c), though this should be confirmed quantitatively. 



Table 1. Effective initial gluon densities, dN 9 /dy, |28I32| and 
time-averaged transport coefficients (q) |35|. for the strongly 
interacting media produced in central A+A collisions at four 
different center-of-mass energies obtained from parton energy 
loss calculations reproducing the observed high pr leading 
7T° suppression at each collision energy ("jet tomography"). 
The measured charged particle multiplicity densities at mid- 
rapidity [45) . dN^ p /drj, are also quoted for each ^/s NN . 





(GeV) 


dN B /dy 


<<?> 

(GeV 2 /fm) 


dN^/dr, 


SPS 


17.3 


400 


3.5 


312 ± 21 


RHIC 


62.4 


650 


7. 


475 ± 33 


RHIC 


130. 


~900 


~ll 


602 ± 28 


RHIC 


200. 


1100 


14. - 15. 


687 ± 37 



5 QCD medium properties via 
tomography" 



"jet 



Fig.Elcompiles all the available Raa(pt) for highpy (lead- 
ing) neutral pions measured in central A+A collisions in 
the range ^/s WJV « 20 - 200 GeV. The experimental sup- 
pression factors can be well reproduced by parton energy 
loss calculations that assume the formation of strongly 
interacting systems with initial gluon densities per unit 
rapidity in the range dN 3 /dy w 400 - 1200 [251152] or, 
equivalently jxj , with time-averaged transport coefficients 
(q) « 3.5 - 15 GeV 2 /fm US] (see TableEJ. 



High p T 7t° in 0-1 0% central A+A: 

^ Pb+Pb @ \Js^, = 1 7.3 GeV (WA98) 

Au+Au @> = 62.4 GeV (PHENIX Preliminary) 
Au+Au @ = 130 GeV (PHENIX) 

Au+Au @ \|si^ = 200 GeV (PHENIX Preliminary) 
GLV energy loss (dN*/dy = 400) 
GLV energy loss (crf/dy = 650) 
GLV energy loss (dN°/dy = 1 100) 
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Fig. 6. Nuclear modification factor, Raa(pt), for highpr pion 
production in central nucleus-nucleus reactions in Pb+Pb )18| 
at v /s^ = 17.3 GeV; and Au+Au at ^/s~^ = 62.4 GeV [IH] . 
v /s^77 = 130 GeV HI], and v /s^7 = 200 GeV |55| : compared 
to GLV parton energy loss calculations [2,8)13. 21 for different 
initial gluon densities (dN 3 /dy = 400, 650 and 1100). Absolute 
normalization errors of the data, O(10%-25%), are not shown. 



For each collision energy, the derived values for the 
initial rapidity density, dN 9 /dy, and transport coefficient 
(q) , are consistent with each other and with the final par- 
ticle density measured in the reactions. Indeed, assuming 
an isentropic expansion process, all the hadrons produced 
at mid-rapidity in a A+A collision come directly from the 



original gluons released in the initial phase of the reac- 



tion 



dN 9 N t , 



dy 



N r . 



dr) 
dy 



dNch 
drj 



1.8- 



dN, 



ch 



dr) 



(4) 



This relation is relatively well fulfilled by the data as can 
be seen by comparing columns third and fifth of Table ^ 
The time-dependent transport coefficient scales with the 
energy density of the medium 5 (e in GeV/fm 3 ) like [36) : 



q{r) m 8 • £ 3/4 (r). 



(5) 



Since, for an ideal QGP (with 2+1 flavors, i.e. degener- 
acy g w 42), the particle (p f» 4.7 • (T/hc) 3 ) and en- 
ergy (e w 14 • T 4 /(fic) 3 ) densities 03] are related via 
p sa 0.66 ■ (e/hc) 3 / 4 ~ 2.3 • e 3 / 4 , one can express Eq. © as 



q{r) « 3.5 • p(r) 



3.5.p ( 7 j=3.5.— 



7) • <«> 



where for the second equality we have assumed a 1-dim. 
Bjorken expansion. In this scenario, since the medium ex- 
pands boost-invariantly in the longitudinal direction, we 
can further write dV = At tq dy where At is the trans- 
verse area of the system, and therefore 



dN 9 
dy 



(7) 



According to the relation between the time-averaged 
q(r) in an expanding medium and that of a fixed static 
medium is (taking Tf 3> tq): 



(«(r)> 



TO+L e ff 

(T-T )q(T)dT 



2 3.5 dN 9 



^-T~ (8) 



L eff A T dy 



where L eJ y is the effective length traversed by the parton 
in the medium. Despite the simplifying assumptions used, 

4 We use here: N to t/N ch = 3/2, and \dy/dy\ w 1.2. 

5 Rather than a thermodynamical variable, q is actually a 
dynamical quantity resulting from the product of the time- 
dependent density of scattering centers times the strength of 
each single elastic scattering 0. 
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this approximate relation between the medium transport 
coefficient and the original gluon rapidity density is rel- 
atively well fulfilled by the data too (Table [I}. E.g. by 
taking L e g w 4 fm and (At) « 125 fm 2 for 0-10% central 
Au+Au we get 



(q(r)) w 0.014 



dN 9 
dy 



(9) 




6 Excitation function of high p r leading 
hadron suppression 

Based on rather general grounds, the total amount of lead- 
ing hadron suppression at a fixed (large) pt in central 
A+A collisions should depend on the collision energy only 
via two -^i-dependent factors: 

(i) the initial parton density of the produced system, and 

(ii) the relative fraction of quarks and gluons fragmenting 
into the hadron at the pt value in question. 

Indeed, on the one hand, since AEi oss oc dN 9 /dy oc dN c h/o 
and since the total particle multiplicity produced at mid- 
rapidity in A+A collisions is observed to follow the ap- 
proximate scaling |45j : 

dN ch /dr) » 0.75 • (N part /2) ■ hx(^/I^/1.5), (10) 

one expects the amount of suppression to increase accord- 
ingly with ^s NN (note, however, that the true evolution of 
the suppression will be faster than that given by Eq. lflT)|) 
since, for increasing energies both dN 9 j dy and the corre- 
sponding lifetime of the quenching medium are larger). 
On the other hand, the probability for a gluon to lose en- 
ergy is a factor Ca/Cf — 9/4 larger than the probability 
for a quark 6 , and the relative fraction of hard scattered 
quarks and gluons (going through the medium and) frag- 
menting into a hadron at a fixed px varies with ^Js NN 
in a proportion given by a tradeoff between (i) the rel- 
ative density of quarks and gluons at the corresponding 
Bjorkcn x = 2pr/\/~s, and (ii) the relative fragmentation 
"hardness" of quarks and gluons at the corresponding z 
value. A full NLO calculation [47] gives the results shown 
in Fig.|3 (bottom). 

In reference fJHj , Wang&Wang presented a pQCD cal- 
culation of the expected ^/s^^-dependence of the nuclear 
modification factor for high-p^ k° production in central 

6 In QCD, the relative strengths of the three distinct quark 
and gluon vertices, a a CF for q — > qg, o s Ca for g — > gg, and 
u s Tf for g — > qq are completely determined by the struc- 
ture of the gauge group (Casimir factors) describing the strong 
force. For SU(N C ) where N c is the number of colors, Ca = N c , 
C F = (N 2 - 1) /2N C and T F = 1/2. The probability for a gluon 
(quark) to radiate a gluon is proportional to the color factor 
Ca = 3 (Cf = 4/3). In the asymptotic limit, and neglecting 
the splitting of gluons to quark-antiquark pairs (proportional 
to the smaller color factor Tr = 1/2), the average number of 
gluons radiated by a gluon is, therefore, a factor Ca/Cf = 9/4 
higher than the number of gluons radiated by a quark |4ri) . 
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parton -* hadron+X at p T = 5 GeV/c : 




Fig. 7. Top (Fig. taken from |48|1: Raa for neutral pions as 
function of collision energy at fixed pr — 6 GeV/c in 10% 
most central AuAu collisions for non-Abelian (lower curve) and 
"non-QCD" (upper curve) energy loss patterns. Bottom: Rela- 
tive proportion of quarks and gluons fragmenting into a hadron 
at fixed pr = 5 GeV/c in p+p collisions in the range i/s — 10 
- 5500 GeV as given by NLO pQCD g7]. 



Au+Au collisions due to parton energy loss in the pro- 
duced (2-D expanding) QGP. The resulting curve is shown 
in Fig. (top). The amount of suppression (for a 6 GeV/c 
leading hadron) increases monotonically with ^s NN due 
to the growing initial parton density, QGP lifetime, and 
gluonic nature of the quenched parton, and seems to sat- 
urate at Raa ~ 0.02 for cm. energies above ~3 TeV. The 
existence of a maximum amount of suppression is due to 
"irreducible" particle production from the outer corona 
of the medium, which remains unsuppressed even for ex- 
treme energy densities |3rJ] . 

In order to test the effect of the radiative QCD energy 
loss, they compared the expected non-Abelian prescrip- 
tion (in which gluons lose AE q /AE g =9/4 times more 
energy than quarks) to an arbitrary "non-QCD" recipe in 
which quarks and gluons lose the same amount of energy 
(AE q — AEg). For a fixed hadron pr value, say pr ~ 
6 GeV/c, the total suppression factors from non-Abelian 
and non-QCD energy losses are relatively similar below 
^/s NJV 100 GeV, since quarks are the dominant par- 
ton fragmenting into a high px hadron (Fig. bottom). 
Above y/s NN w 100 GeV, gluons take over as the dominant 
parent parton of hadrons with px ~ 6 GeV/c and, conse- 
quently, the Raa values drop faster in the canonical non- 
Abelian scenario. The experimental excitation function of 
high pr 7r° suppression in central A+A collisions sup- 
ports the expected QCD radiative energy loss behaviour 
as demonstrated in Fig. |H1 
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Non-Abelian energy loss: A E g /A E q = 9/4 

"Non-QCD" energy loss: A E g = A E q 

n 1 

u 20 40 60 80 100 120 140 160 180 200 220 

\Tsnn (GeV) 

Fig. 8. Excitation function of the nuclear modification fac- 
tor, Raa(^/s nn ), for 7T° production in central A+A reactions 
at a fixed pr = 4 GeV/c value, compared to predictions of 
a jet-quenching model with canonical non-Abelian (solid line) 
and "non-QCD" (dashed line) energy losses |48|. The shaded 
band around each data point represents the absolute system- 
atic errors (absolute normalization of A+A and p+p spectra 
and nuclear overlap, (Taa), uncertainties) 

Summary 

Experimental results on high px leading hadron produc- 
tion in nucleus-nucleus reactions at center-of-mass ener- 
gies y/s NN « 20 - 200 GeV have been discussed with an 
emphasis on the observed suppression of the per-nucleon 
yields relative to p+p collisions at the same y/s. The amount 
of suppression steadily increases from CERN-SPS ener- 
gies (Pb+Pb at y/s NN « 20 GeV) reaching a maximum 
quenching factor of ~5 at the highest RHIC energies (Au+A 
a t \/ s nn = 200 GeV) due to the increased initial parton 
density, lifetime of the dissipative medium, and gluonic na- 
ture of the parent fragmenting parton. The pt and y/s NN 
dependences of the measured nuclear modification factors 
are in agreement with theoretical calculations of final-state 
non-Abelian energy loss in a dense QCD medium. The 
observed dependence of the suppression factors on the re- 
action plane orientation is, however, significantly stronger 
than the one expected from the oc L in-medium path- 
length dependence predicted by jet quenching models, and 
points to an additional source of azimuthal anisotropy 
(likely collective elliptic flow) in hadron production at pt 
values above 4 GeV/c. 
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